Critical Zeeman Splitting of Fermi Superfluidity at Infinite Scattering Length 
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We determine the critical Zeeman energy splitting for Fermi superfluidity at infinite s-wave scat- 
tering length according to the Monte Carlo and experimental results of the equations of state. Based 
on the universality hypothesis, we show that there exist two critical fields H c \ and H C 2, between 
which a superfluid-normal mixed phase is energetically favored, and model-independent formulae 
for H c i, H C 2 and the critical population imbalance P c are derived. Using recent Monte Carlo and 
experimental results of P c , H c \ and H C 2 are determined. It is found H c i = 0.41ep and H C 2 — 0.50eF, 
with £f being the Fermi energy of non-interacting gas. 
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While Bardeen-Cooper-Schiffer (BCS) superconduc- 
tivity/superfluidity in Fermi systems has been investigted 
many years ago[l|, the main scientific interest in recent 
experiments of two-component ultracold Fermi gas is to 
create Fermi superfluids in the BCS-Bose-Einstein con- 
densation (BEC) crossover [2]. At a wide Feshbach res- 
onance point where the s-wave scattering length a s di- 
verges, a novel type of Fermi supcrfluid has been ob- 
served. In the dilute gas limit, there exists only a single 
momentum scale in the problem, namely the Fermi mo- 
mentum fcp = {S^n) 1 ^ 3 which reflects the inter-particle 
distance n -1 / 3 . In this so-called unitary Fermi gas, any 
physical quantity can be expressed in terms of its value 
in the non-interacting case multiplied by a universal con- 
stant For example, the energy density can be written 
as £ = ££o, where £ is a universal constant and £q the 
energy density of non-interacting Fermi gas. 

In addition to the idealized case where fermion pairing 
happens on a uniform Fermi surface, the effect of Zeeman 
energy splitting Ez = HbH between spin- up and -down 
electrons on BCS superconductivity was known many 
years ago (in the following we absorb the magnetic mo- 
ment /xb in to the definition of the "magnetic field" H). 
At a critical Zeeman field or the so-called Chandrasekhar- 
Clogston (CC) limit H c = 0.707A where A is the 
zero temperature gap, a first order phase transition from 
the gapped BCS state to the normal state occurs. Fur- 
ther theoretical studies showed that the inhomogeneous 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state [5(] may 
survive in a narrow window between H c and -Hfflo — 
0.754Ao. However, since the thermodynamic critical field 
is much smaller than the CC limit due to strong orbit 
effect 0|, it is hard to observe the CC limit and the FFLO 
state in ordinary superconductors. 

Recent experiments for strongly interacting ultracold 
Fermi gases give an alternative way to study the pure Zee- 
man effect on Fermi superfluidity [f| . The atom numbers 
of the two lowest hyperfine states of 6 Li atom, denoted 
by JVf and N^, are adjusted to create a population im- 
balance or polarization P — (jVj — N\)/ (JVj + iVj,) , which 
simulates effectively the Zeeman field H in a supercon- 
ductor. At unitary, phase separation between unpolar- 



ized superfluid and polarized normal gas, predicted by 
early theoretical works [^, [l(| and Monte Carlo (MC) 
simulations [ll|, is observed, but the evidence for the 
FFLO and the breached pairing 0, @] states is lacked. 

An important scientific problem is to determine the 
ratio between the critical Zeeman energy splitting and 
the Fermi energy for a homogeneous unitary Fermi gas, 
H c /e-F, which should be a universal constant at unitary. 
In an early work ll|, the CC limit at unitary was pre- 
dicted to be 
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where — A(/i)//i with /i being the fermion chemical 
potential is another universal constant at H = 0. With 
the MC data of £ and /3, H c /A ~ 1 is found. However, 
this result denotes only the first order superfluid-normal 
phase transition point in the grand canonical ensemble 
with the chemical potential /i fixed, and hence not the 
wanted result H c /cf. On the other hand, this result is 
obtained assuming the normal phase is fully polarized, 
but recent MC work [l2J and experiment [l3[ show that 
the normal phase at the phase transition is partially po- 
larized. 



In addition, in Ref. 



14 , [15( it is shown that 
there exist two critical fields H c i and H C 2 in the BCS- 
BEC crossover and phase separation is the energetically 
favored ground state in the region of H c \ < H < H C 2, like 
the type-II superconductors. Since the particle numbers 
ATj and N± are used as tunable parameters in MC calcu- 
lations and experiments, only the first order phase tran- 
sition point (H/n) c and the critical spin population P c 
are directly determined, and the two CC limits H c i,H C 2 
for homogeneous Fermi gas has not yet been measured in 
MC calculations [12| and experiments [13| |. 

In this paper, we will determine the two CC limits H c i 
and H c2 in terms of the Fermi energy ep for a homoge- 
neous Fermi gas at infinite scattering length, based on 
the universal property of the thermodynamics . While 
a similar approach was considered in [Hill, the mam 
purpose of our work is to determine the two critical fields 
for homogeneous Fermi gas, and our study should be ex- 
act, since we do not assume the normal phase is fully 



2 



polarized. 

At unitary, we can construct the exact equation of 
state(EOS) in the grand canonical ensemble from the uni- 
versality hypothesis. The pressure in the polarized nor- 
mal phase (N) as a function of averaged chemical poten- 
tial fi — and the Zeeman field H = (/if— A*j.)/2 
takes the form [3| 
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where fi^ and /it are the effective chemical potentials 
for the two spin components, M is the fermion mass. 
G{x) is an unknown universal scaling function, and only 

— £n 3 ^ 2 i s known to be £n — 0.56. The to- 
tal number density and magnetization are, respectively, 
n^(fj,,H) = dVu/dfJ, and m^(fj,,H) = dV^/dH. There 
exists a universal critical value So of H//i, below which 
the gas is in the partially polarized state (Npp) with 
ton < tin- In the fully polarized normal state (Npp) 
with run = riN, the gas should be non-interacting with 
Q{x) = ^(1 + xf' 2 . While mean field theory predicts 
So = 1 [14| . it is 3.78 11211 from recent MC simulations. 
Monte Carlo studies[13 and experiments 13] show that 



the stable superfluid phase (SF) at unitary should be un- 
polarized. The pressure does not depend on H explicitly, 
and takes the well-known form 



b/2t-3/2 



(3) 



The total density reads ns F (/U, i?) = cfj, 3 ^ 2 ^ 3 ^ 2 . 

In the grand canonnical ensemble, a first order SF-N 
phase transition occurs when the pressure of the super- 
fluid and the normal phase equate, i.e. , when H j [i reaches 
another universal critical value 7, which is determined by 
the algebra equation 



e(7) = r 3/2 . 



(4) 



If the normal phase at the phase transition is Npp, we 
have + H/fi) 5/2 = £~ 3/2 . Combining the relation 
A(/x) = /3/x at unitary, we immediately recover the result 
© derived in [ll|]. 

To determine the critical fields H c \ , H C 2 for homoge- 
neous Fermi gas, we turn to the canonical ensemble with 
fixed total particle number density n = (2Me F ) 3 / 2 / (37r 2 ). 
The Zeeman splitting H is treated as a real external field, 
and the conversion between particles in the states f and 
I is allowed, but the chemical potential fj, is then not a 
free parameter. 

The chemical potential fi^(H) in the polarized normal 
phase is solved from the number equation %(/in, h) = n, 
and the energy density £^(h) = ^N(H)n—V^(^(H), H) 
reads 
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with £q = |ce F ^ 2 . The well-known relation £ — 3V/2 Q 
breaks down here, since the interacting energy with ex- 
ternal field H is included. In the phase Npp, we have 
fj,jss(H) = 2 2 / 3 ep — H. The Npp-Npp transition occurs at 
H = 2 2 / 3 e F 5o/{l + S ). While mean field theory predicts 
H = 2~ 1 / 3 e F ~ 0.794e F , we find H ~ 1.26e F from the 
recent MC simulations. 

Solving the number equation ns F (/i, H) — n for the 
superfluid phase, the chemical potential and energy den- 
sity are given by hsf(H) — £e F and £s F (H) = ££o- At 
H = 0, there is the BCS instability £ SF (0) < £ N (0) 
which is numerically supported by the fact £n > £■ 
While £g F (H) keeps independent of H, £^(H) should 
be a monotonously decreasing function. If there exists 
no heterogeneous mixed phase, a phase transition occurs 
at £n(H c ) — £gp. Once Q(x) is known, one can deter- 
mine H c . If we assume the normal state at H c is fully 
polarized, i.e., H c > Ho, we find 
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Here Ao = /3£cf is the energy gap at fixed density. From 
H c > H , there is £ < 2 2 / 3 /(l + S ) ~ 0.33, which is 
inconsistent with recent MC result £ ~ 0.42 or £ ~ 0.44. 
On the other hand, taking £ ~ 0.4 — 0.5 and (3 ~ 1.2 
from the recent MC data, we have H c /Aq > 1.5, which 
is in contrast with the constraint H c < Ao [161 ] to ensure 
the superfluid phase is unpolarized. Therefore, the MC 
data on £ and [3 indicates that the phase transition is not 
directly into the state Npp. 

The absence of heterogeneous SF-N mixed phase in 
above analysis is not adequate since the first order phase 
transition is often associated with the phase separation 
phenomenon at fixed density. We then take the heteroge- 
neous mixed phase into account. In this case, the critical 
field H = H c (n) = 7/i for the first order phase transition 
in the grand canonical ensemble splits into a lower and 
a upper critical fields H c \ and H C 2, and the SF-N mixed 
phase appears in the region H c \ < H < H C 2- H c ± and 
H C 2 can be determined by equating the chemical poten- 
tial \x to its value in the superfluid and the normal phase 
respectively, H c \ = 7/Usf an d H C 2 = 7^n(#c2), where 
(J>n{H C 2) is obtained by the number equation in the nor- 

mal state, f i N (H c2 ) = £e F [l - ^ 3/2 lG'(l)] ■ Thus 
we find the following model-independent expression for 
the lower and upper critical Zeeman fields 



Hd = J&f, H c2 = 7^£f 
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From the phase equilibrium condition Vs F (^, H) = 
■Pn(M) H)> m the mixed phase the ratio H / ' \i keeps a con- 
stant 7 and the chemical potential reads (ips(H) = H/j. 
The property of the normal bubble in the mixed phase 
depends on the value of 7. The MC calculations [l2| pre- 
dict 7 < So, i.e., the normal bubble is partially polarized. 
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The volume fractions of the superfiuid and normal phase, 
denoted by x and 1 — x respectively, are determined by 
n = x(H)n SF (^i PS , H) + [l-x(H)]n N (fip S , H). Using the 
EOS for the phases SF and N, we find 



x(H) 



(H/H, 
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The energy density of the mixed phase, £ps(H) — [ipsn- 
'Ppsiups, H), can be evaluated as 



H 

H c i 



3/2' 
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Since the normal bubble is polarized, there is a nonzero 
global polarization P in the mixed phase, i.e., the system 
becomes "spontaneously megnetized" when H > H c \. 
From the definition P = (A T - N l )/(N^ + N{), we find 



P{H)= 1 ~ 1 (H/H* 
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The mixed phase continuously link the superfiuid and 
normal state. We have /ips = fJ-SF at i/ = H c \ and 
MPS = Mn at H = H C 2, which ensures < x < 1 with 
x(H c i) = 1 and x(H C 2) = 0. 

There are three constraints on the universal constants 
and scaling function. The appearance of the mixed phase 
requires H c \ < H c2 , which gives rise to < 7^' (7) < 
|£— 3/2 Since the superfiuid phase is unpolarized, there 
should be H c \ < Ao or f3 > 7, which plays the same 
role as the lower bound for the ratio fi-\ / /x j proposed in 
Finally, the normal state at H = -ff C 2 is partially 
polarized, there is H C 2 < Hq. 

For the discussions above, the mixed phase is assumed 
to be the ground state in the region H c \ < H < H C 2- 
As a complete study, we have to prove that the mixed 
phase has the lowest energy in this region. While the 
information on the scaling function Q(x) is still lacking, 
the requirement of the lowest energy can tell us some 
constraints on it. 

Firstly, the energy density of the mixed phase can be 
written as £ PS (H) = ^£ f{H/H cl ) with f(z) = z - 
|z 5 / 2 . Since f'(z) = 1 — z 3 / 2 , £ps(H) is a monotonously 
decreasing function of H in the region H c \ < H < H C 2. 
Combining with the fact that £ps = <?sf at H — H c \ and 
£sf is H- independent, there is always £pg(H) < £sf(H) 
for H c i < H < H c2 . 

Secondly, the condition £ps(H) < £~n(H) requires 
.9(7) < 5(7') with g(t) = h/t - lG(t) (h/tf 2 , h = H/e F 
and 7' = H/fi^(H). Even though the full information of 
the scaling function Q(x) is lacked, it is sufficient to show 
£ps(H) < £-n(H) at H < H C 2, due to the continuity and 
the fact of £n(0) > £sf(0). From the first order deriva- 
tive of g(t) at t = 7, ff '( 7 ) = l~ 2 h [(H/H c2 ) 3 / 2 - l], 
g(t) is a decreasing function near t = 7. Therefore, at 
H < H C 2 the condition 3(7) < 3(7') requires 7' < 7 or 



Mn > MPS- Since /Un(0) > A*sf and hps is an increasing 
function of H, we believe the relation /xsf < MPS < Mn 
holds in the region H c \ < H < H C 2- A schematic plot of 
the energy density for various phases are shown in Fig[T] 



UJ 
UJ 




FIG. 1: The energy for the superfiuid, normal and mixed 
phases in the region H c i < H < H C 2- 



While the universal constant £ and 7 has been de- 
termined in MC calculations and experiments, to deter- 
mine the upper critical field H C 2 we need the value Q' (7), 
which is not known so far. However, another important 
quantity, the critical population imbalance P c — P(H C 2) 
is determined [H, El- We now consider systems with 
separately fixed N-\ and and without conversion be- 
tween the up and down states. From P(H c i) = and 
P(H C 2) = P c , the ground state is the unpolarized super- 
fluid at P = and the SF-N mixed phase for < P < P c . 
In the mixed phase, the effective "magnetic field" is given 
by H(P) — 7£ep(l + 7P) 2 / 3 , and the critical population 
imbalance P r reads 



Pr 



2G'(7) 
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(11) 



It is very interesting that to determine the values of 
H c i,H c2l it is not necessary to know the full imforma- 
tion of G(x), we need only the value of Q'(x) at x = 7. 
On the other hand, we find there exists a simple relation 
between the CC limits in different cases, 



H 



c2 



H c \ 



= (1 + iPc) 
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To show the P c we theoretically obtained above is con- 
sistent with that obtained in MC calculations and exper- 
iments, we derive the energy density £ps as a function of 
the ratio uj 



£ PS (n h n l ) = -n 



(67r 2 n T ) 2 / 3 
2M ' 



(13) 
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With fixed iVj and N±, the energy density defined as 
£ = jLti-n-f + A 1 !"! — V satisfies the relation £ = 3P/2 in 
all phases, since H is now no longer treated as an external 
field. The function I(z) can be shown to be 



/(z) = 2- 2 / 3 £[(l + 7 ) + (l- 7 )z 



5/3 



(14) 



which is consistent with the formula used in the MC 
calculation [l2T| to obtain P c . 

We can now determine the critical Zeeman fields 
H c \, H c2 from the data of £, 7 and P r , and compare them 
with that obtained with mean field [141] and beyond mean 
field theories [Hj]. The MC and experimental data quite 
close to each other. The MC calculation gives £mc = 
0.42(1), 7 mc = 0.967, and P C MC = 0.389, while the ex- 
perimental data are 7exp = 0.95, and P C EXP = 0.36(£ is 
not measured in [HI]). We thus take the MC data to de- 
termine H c i and H c2 . Substituting them into equation 
(fTTj) . we find G'ijMc) = 2.587. With equation ©, we 
obtain 



H cl = 0.407e F , H c 



0.503e F . 



(15) 



Our formulae (|7|) and (jlip are model independent. In 
mean field theory, we have £mf = 0.5906 and Q(x) = 
\ [(1 + xf' 2 Q{\ + x) + (1 - a;) 5 / 2 e(l - a;)] . Numerical 
solution of equation Q leads to 7 mf = 0.8071 and 
£'(Tmf) = 2.9307. Thus we find H^ F = 0.477e F , 
#mf = .693e F and P C MF = O^Vwhich agree well with 
the numerical values obtained in [14j . One finds the mean 
field value of H c2 deviates significantly from our result. 
Also, the result H c2 = 0.693 + 0.087/A + 0(1/N 2 ) ob- 
tained by the large-N expansion method [3] is also not 
consistent with our result. 

In summary, we have presented a model indepen- 
dent calculation of the lower and upper Chandrasekhar- 
Clogston limits of a unitary Fermi superfluid. Future 
studies should focus on the calculation of the scaling func- 
tion Q{x) beyond- mean-field theories [H, [lj| [2(| • Once 
the universal constant £ and function G{x) are known, 
one can directly obtain the critical polarization P c from 
our model independent formula (|11| and check the consis- 
tency between theory and experiment or MC simulation. 

The above model independent approach can be gener- 
alized to finite temperature T, where both the superfluid 
and normal phase are polarized due to thermal excita- 
tions. From the universality, the EOS for the normal 
phase and superfluid read[|| 

7> n .sf(T, h, H) = Icv 5/2 Gn,sf (H/h, T/ M ) , (16) 
5 

where we have set ks = 1- The scaling functions for the 
superfluid and normal phase should be different. 

In the grand canonical ensemble, one expects that the 
phase transition along the T / /i axis is of second order 
at small H/n and first order at large H/fi. The first 
order phase transition is determined by the equation 



Gs{H/n,T/n) = G SF (H/fi,T/fi), or explicitly H/p = 
W(T/(i) with known W(0) = 7. The first order phase 
transition should end at a so-called tricritical point 
(H/ n, T I \x) = (a,b). At mean field level, it is predicted 
to be (a, b) = (0.70,0.38) [Hj]. 

At fixed total particle number, fi is not a free 
variable, and the tricritical point is characterized by 
(Ttcp j Htcp ) ■ Due to the continuity with the zero tem- 
perature case, for T < Ttcp , there exist two critical fields 
H cl (T) = /tiVVfT/zti) and H c2 (T) = n 2 W{T/n 2 ) with 
[L\ and fi 2 being the chemical potentials corresponding to 
H c i and H c2 . The mixed phase region H c i < H < H c2 
should decrease with increasing T, and finally disap- 
pears at the tricritical point with -ffxcp = o2tcp/& & n d 
^1=^2= T TCP /b. 

When JVf and JVj_ are fixed, for T < Ttcp, the phase 
separation should be the ground state in the region Pi < 
P < P 2 with P 1 = P(H cl ) and P 2 = P{H c2 ). At T ^ 0, 
Pi should be nonzero and increase with temperature, and 
Pi = P 2 = Ptcp at the tricritical point. Once the scaling 
function Q and the tricritical point (a, b) are known, Ptcp 
and Ttcp can be calculated from the following model- 
independent formulae: 



Ptcp 

Ptcp 
£f 



g> x (a,b) 



^G(a,b)-aG' x (a,b) 



bg>(a,b) 



5Ptcp 
2<%M 
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(17) 



with the definition Q' x (x, y) = dQ(x, y)/dx and G' v (x, y) — 
dG{x,y)/dy, and G can be the scaling function of either 
the superfluid or the normal phase. 
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